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We derive the mean eigenvalue density for symmetric Gaussian random N x N matrices in the 
limit of large N, with a constraint implying that the row sum of matrix elements should vanish. The 
result is shown to be equivalent to a result found recently for the average density of resonances in 
random impedance networks [Y.V. Fyodorov, J. Phys. A: Math. Gen. 32, 7429 (1999)]. In the case 
of banded matrices, the analytical results are compared with those extracted from the numerical 
solution of Kirchhoff equations for quasi one-dimensional random impedance networks. 



The study of random matrices^ has provided insight 
into many physical problems, both in the quantum and 
in the classical domain. For example, random matri- 
ces have been very successfully used to model statistical 
properties of disordered conductors and of highly excited 
classically chaotic quantum systems 2 . In the classical do- 
main random matrices arise, for instance, in the context 
of diffusion in random, directed environments (see for ex- 
ample Ref. y and references therein) . In most of those 
applications the random matrix elements obeyed symme- 
try requirements where appropriate, and were otherwise 
taken to be independently distributed random variables. 

There are cases, however, where constraints on the ma- 
trix elements must be considered. Such constraints gen- 
erate correlations between the latter: examples are elec- 
tron hopping in amorphous semiconductors (see Ref. 
and references cited therein), random impedance net- 
works^&l (see Ref. for a review and Ref. for ap- 
plications), and random master equations 10 : in these 
cases, the random matrices obey the constraints that 
the row sums of matrix elements should be zero. This 
condition implies correlations between diag onal and off- 
diagonal matrix elements. In Refs. 4 6. lfj, the average 
spectral density and spectral fluctuations of three dif- 
ferent random-matrix ensembles of this type were calcu- 
lated, using the method of replicas and the supersymme- 
try approach. 

In this paper we derive the average density of eigenval- 
ues of a suitably modified ensemble of symmetric Gaus- 
sian random matrices. The aim is twofold. First, we 
wish to show that for full matrices the density obtained 
is equivalent to that found in Refs. lol and Hoi Our second 
aim is to derive the corresponding result for banded ma- 
trices, and to compare it with the results of the numerical 
solution of the Kirchhoff equations on random impedance 
networks with quasi one-dimensional topology. 

Formulation of the problem. We consider an ensemble 



of N x N random matrices M with matrix elements 
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where J is a real symmetric N x N matrix with ran- 
dom entr ies. Such ensembles have been considered in 
Refs.fflnl andlU In Ref.El M was used to model a 
random transition matrix for a model describing glassy 
relaxation, dtu(t) = — M«(t), with matrix elements 
J mn distributed independently (subject to the constraint 

J ran = Jnm) according to 
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The form JQ| yields ^ n M mn = 0, implying probability 
conservation in the problem considered in Ref. Ilfj . The 
average density d(X) of eigenvalues A, 



d(X) = iV- 1 {tr5(M - Al)), 



(3) 



was calculated in the limit of large N and p, using the 
method of replicas. (• • •) is an average over the ensemble 
defined by P(J mn ). 

In Ref. 0, eigenvalues of matrices M with elements 
similar to Q were shown to model resonance frequencies 
in random impedance networks^, with J mn — J nm and 



J* {Jmn) 2 ^\^7nn l) 



1 



(4) 



In Refs. |6| and Q, the average density of resonance fre- 
quencies was calculated, in the limit of large N, using 
a variant of the supersymmetry technique. It was found 
that the result agrees^ with that of Ref. 0, up to a scale 
factor (related to p) and a rigid shift in A [related to the 
fact that J mn > in ©]. 



2 



In the following we calculate the ensemble-averaged 
density of eigenvalues of M, treating both off-diagonal 
and diagonal entries of J as independent, identically 
distributed Gaussian real variables. The corresponding 
symmetric random matrix J belongs to the Gaussian or- 
thogonal ensemble (GOE)i with joint probability density 



P(J)dJ oc exp ( 5— tr J 2 

4cH 



dJ. 



(5) 



The average density of eigenvalues E of such matrices in 
the limit N 3> 1 is given by the semi-circular law d(E) = 
(27R7 2 iV)-V4cr 2 iV- E 2 for \E\ < 2(cr 2 A^) 1 / 2 and zero 
otherwissi^. In the following we ask how imposing a 
" constraint" (that the row sum of matrix elements should 
be zero) modifies the mean eigenvalue density of M with 
respect to that of J. In the limit of large N, the problem 
may be solved using diagrammatic perturbation theory 
(see for instance Ref. y and references cited therein) , as 
shown below. 

Our results may be summarised as follows. In the limit 
of large N, the averaged eigenvalue density for M coin- 
cides with the result derived in Refs. 00- In this limit, 
correlations between diagonal and off-diagonal matrix el- 
ements are found to be irrelevant, and the result can be 
understood in terms of an averaged Pastur equation^. 
Furthermore, we have also considered the case of banded 
symmetric random matrices. This case is of interest for 
random impedance networks with quasi one-dimensional 
topology^ii. Our results are in good agreement with 
those of exact numerical solutions of the Kirchhoff equa- 
tions for such networks. 

Method. The average eigenvalue density may be ob- 
tained from the trace of the averaged resolvent G = 
((Sl-M)- 1 ), 



d(E) = -(jVV) _1 ImtrG. 



(6) 



Expanding the resolvent, Wick's theorem may be em- 
ployed for performing the average, using 

(M tj M mn ) = a 2 (dimSjn + S in Sjm) (7) 

+ SijS mn (N5 im + 1) . 

In the limit of large N, adopting the scaling tr 2 = iV -1 , 
can be simplified to 



(8) 



It is convenient to keep track of the contributions with 
the help of diagrams. To this end, a graphical represen- 
tation of the contraction © is introduced in Fig. [fla). 
The averaged resolvent G turns out to be a diagonal ma- 
trix. Fig. njb) defines a graphical representation for the 
diagonal elements Gjj of G. Terms contributing to Gjj 
are shown in Fig.^c-e). One observes that, in the limit 
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FIG. 1: (a) Contraction (MijMji): the dashed line carries a 
factor of 1/N and the wavy line carries a factor of unity. The 
second term contributes only when i = j. (b) Graphical rep- 
resentation of Gjj = E' 1 and Gjj, the jth diagonal matrix 
element of the average resolvent G. (c) Example of a diagram 
contributing to Gjj . Internal indices are summed over (i and 
k). (d,e) Diagrams of higher order in JV -1 . 
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FIG. 2: (a,b) Self-consistent equations for the average of 
the resolvent, (c) Graphical representation of the auxiliary 
functions Ha. 



of large TV, only diagrams with no intersections between 
dashed or between dashed and wavy lines contribute. 
Thus, to leading order in N, the contribution (c) in Fig.Q] 
must be considered, but not (d) or (e), for example. 

All diagrams contributing to the resolvent to leading 
order in iV _1 may be summed, resulting in a set of two 
coupled equations shown in Fig. j^a) and (b). The di- 
agonal elements Gjj and Hjj are independent of j and 
denoted by G and H in the following. The first equation 
[Fig. Ufa)] contains an infinite sum of diagrams. This 
sum may be performed exactly. It is denoted by g{H~ 1 ), 
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FIG. 3: Density of eigenvalues of exact diagonalisations of 
random matrices of the form (1,5) for a = 1, N = 500, 1000, 
and b = 50, 100 (symbols), together with the prediction I13H 
(solid line). 



where 



J — oo z J 



(9) 



The result for the resolvent scales with a 2 and N as 
^N)- 1 / 2 g{H~ x (cr 2 iV) 1/2 ). Consequently, the equa- 
tions shown in Fig. Eta,b) imply the following self-con- 
sistent equation 

G = {a 2 N)- 1/2 g((E - G) (a 2 N) 1/2 ) , (10) 

equivalent to Eq. (51) in Ref. 6 for the average density 
of resonances in an ensemble of highly connected random 
impedance networks^. 

The form of the simplified contraction JSJ implies an 
interpretation of the result l)10[l in terms of an averaged 
Pastur equation: consider a random matrix M = J + 
V, where J is distributed according to JSJ, and V is a 
diagonal matrix with Gaussian random entries Vk with 
zero mean and unit variance, independent of J mn . For a 
given realisation of V, Pastur's equation^ is G = (El — 
V — G) -1 (see also Ref. and references therein). In 
this equation G is the J-averaged resolvent, keeping V 
fixed. One obtains Eq. (|10|) after observing that G is 
diagonal, by averaging over the matrix elements of V. 
This interpretation implies that, in the limit of large JV, 
the correlations between diagonal and off-diagonal matrix 
elements of M [as seen in Eqs. (|1I7|| ] are irrelevant. 

The above procedure is easily extended to the case 
of banded matrices, also of interest in random impe- 
dance networks^ 1 . In the banded case, (J m nJki) = 
o~ 2 (5 m kO~ni + SmiSnk) [which follows from JSJ] is replaced 
by 



(JmnJu) = cr 2 (\rri - n\) (5 mk S n i + 8 m i5 n k) ■ 
The function a 2 (x) is given by 



a 2 (x) 



g 2 for < x < 6/2, 
otherwise. 



(11) 
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FIG. 4: Topology of a (quasi) one-dimensional lattice with 
periodic boundary conditions. Each site n is connected to its 
b — 4 neighbours. 



The bandwidth of J is thus b. In the limit of large N and 
large b, the spectral density is given by a slight modifi- 
cation of itTUll . 



G = (a 2 b)-^ 2 g((E - G) (a 2 bf' 2 ). 



(13) 



Diagrammatically, the necessary changes are most easily 
derived by letting a 2 = N^ 1 and assuming that b = BN 
(with fixed B 1). Then the wavy line in Fig. QJa) 
acquires a factor of B. Furthermore, the dashed line 
in Eq. © also acquires a factor of B. Hence the self- 
consistency equation in the banded case becomes <|13ll . 
Let us also note that the same formula is actually valid 
not only for b ~ N, but more generally for 1 <C 6 <C N. 

This result implies that densities for different val- 
ues of b can be scaled on one single curve by plotting 
Vbd(E/Vb). In Fig. □ solutions of 1(15)1 for a = 1 are 
compared with results of exact diagonalisations of ran- 
dom matrices with N = 500, 1000 and b = 50, 100. We 
observe a very good agreement. The results confirm that, 
for large N and 6, the average density of eigenvalues is 
independent of N, and that it scales with b as expected. 

In the following we show that 1)13(1 also describes the 
density of resonances for certain random impedance net- 
works with a large, but finite connectivity. 

Random impedance networks. Random networks of 
complex impedances are currently used to model electri- 
cal and optical properties of disordered inhomogeneous 
media™. The most common situation is that of a binary 
composite medium, modeled by attributing a random 
conductance to each bond (x, y) of a lattice, according 
to the binary law: 



(To with probability 
<7i with probability 



P, 

q = l-p. 



(14) 



The homogeneous Kirchhoff equations for the electric po- 
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FIG. 5: Plot of the density of resonances for a quasi one- 
dimensional random impedance network (see text). Data for 
three different ranges (b = 60, 100, and 120) (symbols) are 
compared with the theoretical prediction 1131 . 



tentials, 



can be recast as 

(A Q - \A)V = 0, 
with A = &o /((To — )> and 

(AV) X = £(V y -V x ) 
y(x) 



(15) 



(16) 



y£P(x) 

(A Q V) X = Yl (Vy~V*)> 



(17) 



yeQ(x) 

where y(x) are all the sites connected to site x, whereas 
y G f(x) (resp. y € Q(x)) are those connected by a 
conductance cto (resp. a%), so that A — Ap + Aq, Res- 
onances appear as non-trivial solutions to Eq. for 
< A < 1. 

An efficient algorithm allowing for an exact determi- 
nation of all the resonances of a finite sample has been 
developed in Ref. [5]. We have adapted this algorithm 
to the simplest geometry allowing for long-ranged bonds. 
Each site n of a very long chain is connected to its b neigh- 
bors (n — 6/2, . . . , n — 1, n + 1, . . . , n + b/2), as shown in 
Fig. 01 Periodic boundary conditions are assumed. For 
deffiiiteness we choose p — q = 1/2, so that all the res- 
onances are expected to be located at A = 1/2 in the 
b — ► oo limit. 

Our numerical results are shown in Fig. [S] for very 
long periodic chains with ranges b = 60, 100, and 120. 
For each value of 6, we have accumulated a number of 
resonances of order 10 7 . After rescaling the resonances 
according to A = (1 + E) /2 the density is given by ljr3|) 
with a = 1. A very satisfactory quantitative agreement 
with the theoretical prediction is observed. 
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